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We discuss the results of a search for anomaly free Abelian %n discrete symmetries that lead to 
automatic R-parity conservation and prevents dangerous higher-dimensional proton decay operators 
in simple extensions of the minimal supersymmetric extension of the standard model (MSSM) 
■ based on the left-right symmetric group, the Pati-Salam group and SO(IO). We require that the 

' superpotential for the models have enough structures to be able to give correct symmetry breaking 

' to MSSM and potentially realistic fermion masses. We find viable models in each of the extensions 

and for all the cases, anomaly freedom of the discrete symmetry restricts the number of generations. 



I. INTRODUCTION 



o 
O 

o\ 

Supersymmetry (SUSY) is widely believed to be one of the key ingredients of physics beyond the standard model 
(-H ' (SM) for various reasons: (i) stability of the Higgs mass and hence the weak scale; (ii) possibility of a supersymmetric 
dark matter; (iii) gauge coupling unification, suggesting that there is a grand unified theory (GUT) governing the 
I ' nature of all forces and matter. The fact that the seesaw mechanism for understanding small neutrino masses also 

requires a new scale close to the GUT scale adds another powerful reason to believe in this general picture. 
(-H There are however many downsides to SUSY. For instance, while the SM guarantees nucleon stability, in its 

I— I, supersymmetric version, there appear two new kinds of problems: (i) There are renormalizable R-parity breaking 
operators allowed by supersymmetry and standard model gauge invariance, e.g. 



= Xrjk U L, el + A^fc L, dl + A^^^ u'i ^ . (1) 



^ ' M/^ = X. 4- \'. .. ^ \" 

O 

' Here, Q, L, u'^, (F and e'^ denote the left-chiral quark-doublet, lepton doublet, u-type, d-type and electron superfields. 
I Combination of the last two terms leads to rapid proton decay and present limits on nucleon stability imply (for 
' squark masses of TeV; cf. e.g. 0,0,01): 

Ki.Ki, < 10-'^ (2) 

O ■ 

These terms also eliminate the possibility of any SUSY particle being the dark matter of the Universe. Usually 

_^ ^ assumptions such as either R-parity or matter parity (cf. are invoked to forbid these couplings and rescue the 

, proton (as well as dark matter) stability, (ii) A second, more vexing, problem is that even after imposing R-parity, 
• one may have dimension five R-parity conserving operators such as Kijki Qi Qj Qk Li/Mp. Such operators also lead 
5^ I to rapid proton decay. In fact, present nucleon stability limits put an upper limit on K1121, K1122 ^ 10~^ [3]. In what 
follows we will refer to these operators as Q Q Q L. 

To understand that these problems can be attributed to SUSY, recall that the R-parity breaking terms are forbidden 
in the case of SM by Lorentz invariance in the standard model and the QQQ L is suppressed by two powers of Planck 
mass and hence not problematic. 

In this study, we focus on the problem of baryon number non-conservation and require the model to satisfy the 
following constraints: (i) R-parity symmetry is exact so that it prevents catastrophic proton decay; (ii) R-parity 
conserving dimension five proton decay operators of type QQQ L are either forbidden or suppressed to the desired 
level and (iii) the superpotential of the theory has enough structure for ensuring proper symmetry breaking down to 
the standard model and give fermion masses. We look for symmetries that allow all desired terms in the superpotential 
while forbidding the unwanted terms discussed and satisfying the anomaly constraints. 

The anomaly constraints can be related to triangle graphs involving the not only discrete symmetry and the gauge 
symmetries but also gravity. The reason for requiring anomaly freedom is the following: there is a general belief that 



^ An exception to this statement is the case where the lightest SUSY particle is the gravitino and plays the role of the dark matter. 
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non-perturbative gravitational effects such as black holes and worm holes break global symmetries of nature [5|. In 
an effective field theory language, they are parametrized by Planck scale suppressed higher dimensional operators. In 
fact, the above-mentioned QQQ L type operators are a manifestation of the fact that global baryon number symmetry 
is broken by such effects. On the other hand no hair theorem of general relativity says that any non-perturbative 
gravitational effect must respect the gauge symmetries. Therefore if there is a gauged discrete symmetry in the theory 
that prevents the undesirable terms under discussion, they will be absent even after all non-perturbative effects are 
taken into account. 

How does one ensure that a discrete symmetry is a gauge symmetry ? This problem has been extensively studied 
in the literature in the context of the MSSM [l|, 0, H, U H, H, and the general procedure is to calculate the 
various anomaly equations involving the discrete group with gravity and the gauge group i.e. vanishing oi D g g, D 
anomalies where D stands for the discrete symmetry group in question, g is the graviton and G is the (continuous) 
gauge symmetry on which the theory is based. 

In the context of the MSSM, a discrete symmetry has been identified, dubbed 'proton-hexality' in |2|], that 
contains R-parity as subgroup and forbids QQQ L Q • Remarkably, anomaly freedom of this requires the 
number of generations to be 3 Moreover, it has been shown that (given 3 generations) this is the only anomaly free 
symmetry that allows the MSSM Yukawa couplings and neutrino masses while forbidding the dangerous operators 0]. 
Such symmetries can be extended such as to also forbid the fi term [11]. (For an approach to ensure proton stability 
by flavor symmetries see e.g. [HI.) On the other hand, the charge assignment is different for different standard 
model representations. This raises the question whether nucleon stability can be ensured by discrete symmetries in 
(unified) theories where standard model representations get combined in larger multiplets, and the charge assignment 
is hence restricted more strongly. We therefore seek discrete symmetries ensuring sufficient proton stability in three 
gauge extensions of the supersymmetric standard model: (i) the left-right symmetric model [Tsj , (ii) the Pati-Salam 
model [ij and (iii) SO(IO). All these models incorporate the B — L gauge group, which is generally used in the 
discussion of the seesaw mechanism for neutrino masse s fisl ITgI . [TtI . [isl . IToj (for a review see e.g. [23|) and also 
provides one way to guarantee R-parity conservation [2l|, I22ll23|. Due to the higher gauge symmetry, which must be 
broken down to the SM gauge symmetry, specific terms must be present in the superpotential. This poses constraints 
on the discrete symmetry. 

One main result of the study is a connection between the order of the discrete group and the number of generations 
in all the cases. We give examples of viable models for all the different gauge groups. In our discussion, we follow a 
certain 'route of unification'. We start with the left-right symmetric extension of the supersymmetric standard model, 
proceed via the Pati-Salam model to SO(IO) GUTs and finally comment on how our results might be used in string 
compactifications. 

This paper is organized as follows: in section|TTl we discuss anomaly-free gauge symmetries ensuring proton stability 
in left-right models; we proceed to the Pati-Salam model in section Hill and in section HVl we discuss SO(IO) models. 
We give our conclusions in section IVl 



II. LEFT-RIGHT MODEL AND DISCRETE SYMMETRIES 



In this section, we discuss the left-right symmetric extension of MSSM, i.e. the gauge group is SU(3)c x SU(2)l x 
SU(2)r X \J{1)b-l- This amounts to breaking B — L symmetry of the model by either a, B — L = 2 triplet or a 
B — L = 1 doublet. In Table U (a), we show the assignments of the quarks and leptons and Higgs bosons. Tables U 
(b) and (c) list the Higgs sectors of the doublet and triplet models, respectively. 



(a) MSSM part. 



Field 


quantum numbers 


Q 

L 
L" 
$ 


(3,2,1,+i) 
(3,1,2,-1) 
(1,2,1,-1) 
(1,1,2,+1) 
(1,2,2,0) 



(b) Doublet model. 



Field 


quantum numbers 




(1,1,2,+1) 


r 


(1,1,2,-1) 


X 


(1,2,1,-1) 


X 


(1,2,1,+1) 


s 


(1,1,1,0) 



(c) 


Triplet model. 


Field 


quantum numbers 


A" 


(1,1,3,-2) 


A^ 


(1,1,3,-1-2) 


A 


(1,3,1, +2) 


A 


(1,3,1,-2) 



TABLE I: Assignment of the fermion and Higgs fields to the representation of the left-right symmetry group 
SU(3)c X SU(2)l X SU(2)r x U(1)s-l. (a) shows the MSSM sector, (b) shows the Higgs sector of the doublet model and (c) 

the Higgs sector of the triplet model. 
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A. Left-right symmetric models — doublet Higgs case 

The (wanted) superpotential is given by 

On the other hand, the foUowing couphngs must be forbidden (we suppress coefficients): 

t^unwantcd = L + Q^'' + X + Q'"" + Lx + L'' + L Q + Q X L'^ 

+ L' + X ^ + + i * + X • 



(3) 



(4) 



To study the anomaly constraints in this model for an arbitrary Z^r group, we start by giving the charge assignments 
under Zat to the various superfields (denoted as for the field F in the equation below) and writing down the anomaly 
constraints. The anomaly constraints for 'Zn gg, TI^n [SU(3)c]^, Zat [SU(2)l]^ and TLj^ [SU(2)fi]^ respectively are^ 



^9 [6 (gg + 9q0 + 2 (gi + qL^)\ + 4(7$ + 2 (g^ + (/^e + ^ + ^75^^) = mod N' 

Ng [2{qQ + qQ.)] =0 mod TV, 

^3 [3 qQ + qL]+2q<s> + q^ + qY = mod N , 

Ng [3 ggo + qi^c] + 2 q^ + q^a + q^c = mod N , 

where Ng denotes the number of generations. In the first equation 



(5a) 
(5b) 
(5c) 
(5d) 



TV' = 



N , iV odd , 
iV/2 , N even 



(6) 



which follows from equation (10) of fi\. 

The assignments must be consistent with the superpotential ([3]) and has to forbid the terms in ([4]). We scanned over 
possible non-anomalous ZAr<i2 symmetries, keeping the number of generations Ng as a free parameter. Remarkably, 
the smallest viable Ng we found is 3, and the smallest N that works with 3 generations is 6. An example is shown in 
Table mi 





QQ 






q<s> 


Qx 




Ix 




qs 


1 


1 


5 


5 








2 





4 






TABLE II: A viable charge assignment in the left-right model with doublets. 



By giving vacuum expectation values (vevs) to the fields x'^ s-nd x'^, the symmetry is broken to a symmetry 
under which matter is odd while the MSSM Higgs are even. That is, we have obtained an effective R-parity which, 
although there is a gauged B ~ L symmetry, originates from an 'external' Kg. 



B. Left-right models — triplet case 



The transformation properties of the fields under the gauge group are shown in Table U (a) and (c). In this case, 
the right handed neutrino masses arise from the renormalizable couplings in the theory. We have to forbid L and 
{Q'^y^ L'^. There are many anomaly-free discrete symmetries which do the job. The interesting point is that in this 
case, the minimum number of generations is Ng = 2 with N = 2 for the discrete group. The smallest symmetry that 
works for 3 generations is K3. 



^ We ignore cubic constraints (cf. ^8]). 



Qj Lk Qj Lk 

(a) (b) 



FIG. 1: Effective proton decay operators, (a) shows the usual triplet exchange diagram [24l.l2a|. (b) illustrates that the 
amplitude vanishes if the mass partner of the Higgs triplet does not couple to SM particles [27l |. 

III. PATI-SALAM MODEL 

We now proceed to the Pati-Salam model, i.e. the gauge group is Gps = SU(4)c x SU(2)l x SU(2)fi. The new 
feature of this model compared to the left right model just discussed is that the quarks and leptons belong to the 
same representation (see Table HIT)). 







Field 


quantum numbers 


Field 


quantum numbers 






(4,1,2) 




(4,2,1) 




r 


(4,1,2) 


r 


(4,1,2) 




X 


(4,2,1) 




(1,2,2) 




X 


(4,2,1) 






s 


(1,1,1) 



Field 


quantum numbers 




(10,1,3) 


A= 


(10,1,3) 


A 


(10,3,1) 


A 


(10,3,1) 



TABLE III: Assignment of the fermion and Higgs fields to the representation of the left-right symmetry group 
SU(4)c X SU(2)l X SU(2)r: (a) shows the MSSM sector, (b) shows the Higgs sector of the doublet model and (c) the Higgs 

sector of the triplet model. 



The rest of the discussion is completely analogous to the left-right case and in fact the solutions displayed in Table 
II for the discrete charges apply to the doublet case, and as in the triplet version of the left-right model there are 
many solutions. Note that the triplet version of Pati-Salam model does not have proton decay due to SU(3)c C SU(4) 
invariance. 



IV. SO(IO) GUT AND DISCRETE SYMMETRIES 



We now turn our attention to the discussion of discrete symmetries in SO(IO) GUTs. In SO(IO) models, the dimen- 
sion 5 proton decay operators QQQL have two sources: (i}the higher-dimensional coupling [IGm]^ and (ii) effective 
operators emerging from integrating out Higgs triplets [HjllHl (see fi gur e [J). Triplet masses of the order Mqut appear 
to be too small to be consistent with the observed proton life time '26] . As discussed, the coefficient of the QQQL 
operators have to be strongly suppressed. This requires an explanation of why both contributions to these operators 
are simultaneously small. As we shall see, such an explanation might arise from a simple discrete symmetry. 

The proton decay via Higgs triplet exchange can be forbidden by eliminating the mass term for the 10 multiplet 
{H). However one needs to make the color triplets in H heavy so that coupling unification is maintained. This can be 
done by introducing a second 10-plet [H') such that it forms a mass term with H but the color triplet field in it does 

Therefore, from now on we will consider SO(IO) models with 2 10-plets. 
(i) where B — L \s broken by 16- Higgs fields (cf. Table ITVl (b)) and (ii) where 



27| 



not couple to standard model matter 

We will discuss two classes of models: 
B — L is broken by 126-Higgs fields [28|, |29| (cf. Table |TV] (c)). We consider Abelian discrete (Zat) symmetries and 
require that higher dimensional R-parity conserving leading order AB ^ operators, which in this case are of type 
16^ (where the subscript m stands for matter), are forbidden as are all R-parity breaking terms while allowing all 
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terms in the superpotential that are needed to break the GUT symmetry down to the MSSM. Our focus is on proton 
stabiUty, and we leave other issues such as fermion masses and doublet-triplet splitting for future studies. 



(a) MSSM part. 



(b) 16-Higgs model. 



Field 


quantum numbers 




Field 


quantum numbers 




16i 






16_2 


H 


10_2 




■>Ph 


T62 


H' 


IO2 




A 


45o 






S 


54o 



(c) 126-Higgs 



model. 


Field 


quantum numbers 


A 


1262 


A 


126 2 


$ 


2IO0 



TABLE IV: SO(IO) model, (a) shows the MSSM sector, (b) the Higgs content of the 16-Higgs model and (c) the Higgs 
content of the 126-Higgs model. charges (see text) appear as subscript. 



A. 16-Higgs models 

In this class of models, one has an independent motivation for introducing a second 10-plet coming from doublet- 
triplet splitting. Further, apart from a pair of 16 © 16-Higgses, 45- and 54-plets are required to ensure proper 
GUT symmetry breaking down to MSSM (cf. |30|, l3l|, [3^, |33| ) . The superpotential terms that must be allowed are: 
iprni^mH, {tprn "(p h)'^ / : "^h^Phj i ^^'^ ^^'^ S , where tjjm,H are matter and Higgs 16-plets, H, A, S are 10-, 
45-, 54-plets, respectively (see tables |IV] (a) and (b)). This leads to the following constraints on the charges: 



+qH = mod N , 2q^^ + 2q^^ = mod N , (7a) 

I^Ph + 9^j^ = mod N , qn + Qh' = niod N , (7b) 

2qA = mod iV , 2qs ^ mod N , (7c) 

3qs = mod N , 2qA + qs = mod N . (7d) 

Here, we denote the 'En charge for a field F by qp. Next, we list the anomaly constraints, 

W{Ngq^^+q^^ +q:;i;^) + 10 iqH + qH')+ 4:5 qA + 5Aqs = mod iV' (8a) 

2Ngq^^+2q^^ +2q:;i;^+qH +qH' +8qA + l2qs ^ Q mod iV . (8b) 

To forbid the dangerous couplings V'm ipm, H', V'™, ipm ''Ph and V'm ipm ipH H, ipm ipH H' and ipmipn^^ the values 
of the Zjv charges have to be chosen such that they satisfy the inequalities 

2 +qH' ^ mod N , 4 q^^^^ ^ mod N , (9a) 

li'^ + ^ ^ > 3 + q^jj ^ mod N , (9b) 

q-ip^ + qH'.H + qt/jH 7^ mod iV , 9iA,„ + to + qA mod . (9c) 



The smallest symmetry allowing to fulfill all criteria is and requires Ng — 3. A possible charge assignments is 

QiAm — Ij q^'H — ~2, (7^^ — +2, qn = —2, qH' = +2, (745 = (cf. tables HVl (a) and (c)). This charge assignment allows 

for seesaw couplings and the possibility of fermion masses from couplings of type ipm ipm H. The allowed operator 
2 

4'm V'm i^H contributes to both the fermion masses as well as to the seesaw. We note that the charge assignment is 
such that we have 

3 X generation + vector-like (10) 

under SO(IO) x Zg. The model also eliminates the dangerous proton decay operator QQQ L or operator of type 
/M-p. Wc note that, although there are higher-dimensional gauge (and Ze) invariant operators, they do not 
give rise to proton decay operators for the following reasons: 

(i) any possible combination of ipjj and Tp^ that breaks B — L cannot multiply the B — L neutral proton decay 
operators QQQ L {or u'^ e'^). 
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field 


<l<Prr. 


H 


H' 


A 


A 




1 


1 


-1 


-1 


1 



TABLE V: Z3 charge assignment. 



(ii) any B — L neutral combination is also Kg invariant because it has to involve as many tpn as ip^ fields, as only 
the SM singlets with B — L charge ±1 attain a vev. Therefore the product of the non- invariant combination 
Q Q Q L with the B — L neutral combination of ipn and ip^ fields cannot be invariant. 

Consider, for instance, the operators ]gs-[V'm]^ [V'h]^ and ■j^[ipm]'^ [tpH i^n]- S^'^st operator, the expectation 

value of [tpH]'^ vanishes because of the first argument while the second operator is not Zg invariant. 

We also observe that operators like ipj^ H^, which would lead to small diagonal entries in the H — H' mass matrix, 
do not exist for the same reason. That is, in this model, proton decay is only due to dimension 6 operators. 

We refrain from spelling out the detailed phenomenological analysis of this model. However, our preliminary studies 
seem to indicate that one can achieve doublet-triplet splitting and realistic fermion masses while avoiding proton decay 
problem by extending the Higgs content. A complete analysis of these issues are defered to a future publication. 

Let us also comment that, like in the left-right model with doublets, the Zg symmetry gets broken by the ipn and 
ipfj vevs down to a which forbids W^, i.e. acts as an R-parity. This means that R-parity in this SO(IO) model 
does not originate from B — L. 



B. 126-Higgs models 

We now discuss models where the 16h © 16h get replaced by 126 © 126 - the motivation bein g th at R-parity 
becomes an automatic symmetry. Such models have been extensively discussed in the literature [s^ l35l. [sgI Ist} . In 
our context it means that the last two of the four inequalities in Eq. ([9]) do not exist (see Eq. ([13]) below). Instead 
we have the following set of constraints on the charges from anomaly freedom 

m Ng q^^ + 10 {qH + qw) + 126 {qa + g^) ^ mod iV' , (11a) 
2Ngq^^^+'35{qA + qs) +qH + qH' = mod TV . (lib) 

The superpotential constraints can be decomposed in analogs of Eq. ([7|) 

2q^^ +qH = mod N , (12a) 

2gv™ + '?A = mod N , (12b) 

QA + qs = mod N , (12c) 

and analogs of Eq. 

2qii^+qH' + modAf, (13a) 

4(7,0,„ / mod N . (13b) 

Typically in a class of these models, there are only 210 dimensional representations that need to couple to A and 
A fields among themselves as well as with 10 Higgs [H, [s^. These imply that the discrete charge of 210 vanishes 
and also that of A and A are opposite. Substituting these conditions into Eq. (fTTjl . it becomes clear that if there is 
only a single 10 Higgs in the model (or, if qi^i — 0), the requirement of anomaly freedom becomes rather constraining. 
We could only find a Zg symmetry with one generation. However, once we allow for non-trivial qn' , it is possible to 
have simple anomaly free Wijy symmetries which satisfy all constraints, the simplest example being a Z3 symmetry 
with the charge assignment listed in table [V] The Zjv<i2 symmetries with non-trivial qn' require Ng = 3 or larger. 

Another possible symmetry is Sg with the charge assignment listed in table liyi (c) . 

This model has the same effective structure as the 126 model of Ref. [s^, l35l. \3Ql. as far as the discussion of 
fermion masses go (even though it has 2 10-plets, one of the lO's does not couple to fermions due to Zg charge 
assignments). 

For the same reasons as in the 16-Higgs models, the QQQ L type induced proton decay operators are forbidden. 

Again we refrain from a detailed phenomenological analysis of the way this model leads to MSSM at low energies. 
We close the discussion with a comment on how easily the MSSM Higgs fields emerge from the superpotential: the 
relevant part of the superpotential has the form: 

W = Mh HH' + Mq a (A -I- a $i?) + AU + A' $3 + \"A<1>H' + A'" A $ A , (14) 
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which has the right hnear combination of MSSM doublets to maintain all the simple form for the fermion mass results 
of Ref. [33, H^, [la, [13 ■ A detailed analysis of these issues will be given elsewhere. 

C. SO(IO) GUTs in higher dimensions 

Let us now comment on implications of our findings for higher-dimensional models of grand unification, such as 
'orbifold GUTs' [13, [HI [12, [S, [i3, [H, [ill , which provide a simple solution to the doublet-triplet sphtting problem. 
In such models the dimension 5 proton decay can be naturally suppressed [4^ . [43j (while dimension 6 proton decay 
is slightly enhanced [47| ) since here the mass partner of the Higgs triplet has vanishing couplings to matter, as in the 
discussion above. However, (brane) couplings like V'm, ^^^o leading to proton decay, have not been discussed in this 
scheme. A reliable discussion of such operators seems hardly possible in the effective higher-dimensional field theory 
framework. 

One possible way to address this question is thus to embed the model into string theory or, in other words, to 
derive orbifold GUT models from string theory. The first steps for doing so have been performed in Refs. [4R[49l[50|. 
This has further lead to the scheme of 'local grand unificatio n' Isil. [53. [ssl 54 1 . which facilitates the construction 



of supersymmetric standard models from the heterotic string [52l. l53l. l55l |56||. Here, the two light MSSM matter 
generations originate from 16-plets localized at points with SO(IO) gauge symmetry. Some of these models can have 
an R-parity arising as Z2 subgroup of a gauged, non-anomalous B — L symmetry [ssl. [56l. [STj . like in ordinary GUTs 
(however, without the need for 126-plets). On the other hand, QQQL operators remain a challenge [5l,[5g|. The 
fact that these operators could be eliminated so easily in conventional GUTs by simple symmetries leads to the 
expectation that similar symmetries will be helpful in the string-derived supersymmetric standard models with (local) 
GUT structures. One lesson which one might learn from our analysis is that one may derive an effective R-parity and 
suppress QQQL by a discrete (possibly Zg) symmetry under which matter 16-plets have a universal charge. One 
might further hope to get insights about the origin of the discrete symmetries (which remains somewhat obscure in 
the 4D field-theoretic approach) in string models. These issues will be studied elsewhere. 

V. CONCLUSION AND COMMENTS 

Motivated by the beauty of the ideas of supersymmetry and unification, wc have started a search for discrete 
symmetries that forbid proton decay operators in gauge extensions of the supersymmetric standard model. We required 
the symmetries to allow the standard interactions and to be anomaly free. Considering the left-right symmetric, Pati- 
Salam and SO(IO) GUT models with various Higgs contents, we could identify (surprisingly simple) symmetries 
that satisfy all our criteria. In many cases, there is a connection between the anomaly freedom and the number of 
generations. Often, simple symmetries exist only for 3 generations (or multiples thereof), as in [l]. In the SO(IO) 
models, our symmetries forbid dimension 5 proton decay operators. 

Our findings can be interpreted in the following way. Supersymmetric models with an extended or GUT symmetry 
are often challenged by proton decay. To rectify this, one might be forced to introduce additional (discrete) symmetries. 
Our examples then show that R-parity can be a consequence of these additional symmetries rather than being related 
to B — L. From this one might conclude that the appearance of fields with even B — L charges is not a necessity, and, 
for instance, 16-Higgs and 126-Higgs SO(IO) models can be on the same footing. 

It is also interesting that the minimal 126-based SO(IO) models become free of all dangerous proton decay operators 
without losing their ability to be predictive in the fermion sector once we add a simple anomaly free discrete symmetry. 
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